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' ^ '1 Abstract. The dynamics of loops at the DNA denaturation transition is studied. A 

CJ ' scaling argument is used to evaluate the asymptotic behavior of the autocorrelation 

^ . function of the state of complementary bases (either open or closed). The long-time 

G '. asymptotic behavior of the autocorrelation function is expressed in terms of the entropy 

exponent, c, of a loop. The validity of the scaling argument is tested using a microscopic 
model of an isolated loop and a toy model of interacting loops. This suggests a 
method for measuring the entropy exponent using single-molecule experiments such 
C^ . as florescence correlation spectroscopy. 
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1. Introduction 

Melting, or tliermal denaturation of DNA, is the process by which the two stands of the 
DNA molecule become fully separated upon an increase of the temperature [H 12] . At 
low temperatures the strands are partially unbound by forming fluctuating loops where 
the two strands are locally separated. As the melting temperature Tm is approached 
the average loops size increases, yielding full denaturation at Tm- Melting of DNA 
has been extensively studied over the years both theoretically and experimentally. The 
natural order-parameter of the denaturation transition is the fraction of bound base- 
pairs. This was measured using specific-heat and UV absorption experiments. Simple 
models have yielded theoretical expressions for thermodynamic properties. Two main 
approaches have been developed. One, known as the Peyrard- Bishop model, considers 
the two strands as directed polymers interacting via a short-ranged potential [3]. One 
then focuses on the distance between complementary pairs as the melting transition is 
approached. The other, known as the Poland-Scheraga (PS) model, represents the DNA 
molecule as an alternating sequence of bound segments and open loops, and focuses on 
the fraction of bound base-pairs [4]. Within the PS approach self-avoiding interactions, 
which are inherently long-range, may be taken into account. As have been shown these 
interactions affect the loop entropy, which controls the nature of the melting transition 

[HIE]. 

Recently, single-molecule techniques such as optical tweezers [H [8] , magnetic traps 
[HI [in] and Fluorescence Correlation Spectroscopy (FCS) [HI [12], [13] have been used 
to probe properties of the melting process. Other techniques, such as quenching, 
have also been applied [H]. Some of the experiments utilize an external force to 
induce unzipping of the two-strands and study their dynamics. In others, the distance 
between two complementary base-pairs is probed by FCS without applying an external 
force. These experimental methods enable one to study not only bulk properties but 
rather microscopically fluctuating quantities. Inspired by these experiments, theoretical 
treatments of dynamical properties of DNA have been developed. Several studies have 
focused on the dynamics of isolated loops away from the melting transition [151 [IS] ^^^ 
at the transition [17] . The survival probability of an isolated loop has been calculated. 
A toy model for the dynamics of interacting loops has also been introduced and analyzed 

It has recently been shown that studying the loop dynamics may yield information 
on the loop entropy [17]. Within the PS approach the dependence of the entropy of a 
loop on its length plays a dominant role in determining the thermodynamic behavior 
near the transition. On general grounds one can argue that the entropy of a loop of 
length n takes the form 5* = kB\og{Vt{n))^ where Vt{n) ~ s^/rf is the number of loop 
configurations. Here s is a mo del- dependent constant and c is a universal exponent. 
The numerical value of c has been debated over the years. It was found to be modified 
when the excluded-volume interactions, which are long ranged in nature, are taken into 
account [H [5l [6] . When interactions between loops are neglected, and excluded volume 
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interactions are taken into account only within each loop an exponent c ^ 1.76 was 
found [Sj. On the other hand, when excluded volume interactions both within a loop 
and between the loop and the rest of the chain are taken into account, the entropy 
exponent was found to increase to c ~ 2.12 [2l[6]. This latter result, which predicts a 
first order denaturation transition for homopolymers, has been verified numerically |20j . 
While numerical studies of the homopolymer model with excluded volume interactions 
yield a clear first order transition [21] , a direct experimental measurement of c is rather 
difficult and has not been carried out so far. In [17] it was shown that at the melting 
transition the time dependence of the base-pair autocorrelation function depends on the 
parameter c. The base-pair autocorrelation function is defined as Cj(t) = (ui(t+r)Mj(r)) 
where Uj(t) = 1, is a variable which indicates if base pair i is open (1) or closed (0) at 
time t, and (■) denotes an average over r. The behavior of the autocorrelation function 
was studied theoretically away from the melting transition [151 [E] and may, in principle, 
be obtained experimentally by FCS studies. In these experiments the states of a specific 
base-pair is monitored. So far, FCS experiments have been restricted to short molecules 
[TT] . Measuring the exponent c requires extending these studies to longer molecules. 

In the present paper we elaborate on and extend the analysis presented in [1^ for 
the dynamical behavior of homopolymers at the melting transition. The dynamics of a 
single loop is studied using a simple model, whose validity is then verified in detail using 
numerical simulations. Within this model the entropy exponent c is introduced as a free 
parameter which may be chosen at will. While the studies in [TTj were tested numerically 
only for the case c = 3/2, here we test the robustness of the results for models with 
arbitrary values of c. We then consider a toy model, similar to the one considered in 
[TS] . which indicated that the results still hold when the interaction between loops is 
taken into account. 

The paper is organized as follows: In Sec. 2 we study the single loop model using 
both the scaling argument and microscopic models. In Sec. 3 results for the many loops 
model are presented. Finally, we end with a brief summary. 

2. Single Loop Dynamics 

We start by considering the dynamics of an isolated loop. In this approach one ignores 
processes like merging of loops and the splitting of a large loop into two or more 
smaller ones. This may be justified by the fact that the cooperativity parameter, which 
controls the statistical weight of opening a new loop, is estimated to be rather small, 
ctq ~ 10~^ [T9]. Thus splitting a loop into two is unfavorable. Also, the average distance 
between loops, which within the PS model is proportional to l/o"o, is large, making 
the independent loop approximation plausible. In Sec. [3] we introduce a simple model 
to effectively take into account the interactions between loops and show that these 
interactions do not modify the results obtained within the single loop approach. 

Within the single loop dynamics, we assume that a loop may change its length by 
closing or opening of base pairs at its two ends. It survives as long as its two ends do 
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not meet. Let G{nQ, t) be the survival probability of a loop of initial length Uq at time t. 
As discussed above, the quantity of interest is the equilibrium autocorrelation function 

rm ^ T.Z=iPe,{no)n^G{no,t) 

'-'V') v=^3o 5—7 ^ 1 \^) 

E„o=i Peq{no)no 
where for simplicity of notation we have dropped the site index i. Here Peq{nQ) is 
the probability of having a loop of length uq in equilibrium. Hence, rio-Pegl'^o) is the 
probability of a particular site to belong to a loop of length hq. Note that we assume 
that site i remains open as long as the loop survives. This approximation does not affect 
the behavior of the autocorrelation function in the scaling limit. 

We proceed by first presenting a scaling analysis demonstrating that in the case 
of a homopolymer and at criticality, the autocorrelation function decays at large t as 
C{t) ~ t^~'^/^ for c > 2, while it remains finite, C{t) = 1, for c < 2. These results are 
then tested and verified using numerical simulations for various values of c. 

2.1. Scaling Analysis 

In the case of a homopolymer and at criticality it has been shown that the equilibrium 
loop size distribution is Peq{n) ~ 1/n'^. To estimate the survival probability of a loop of 
length no, we consider dynamics under which the loops are non-interacting and do not 
split into a number of smaller loops. Similar to [151 [IHl we further assume that the loop 
is in a local thermal equilibrium at any given time during its evolution. The validity 
of this assumption will be discussed in detail below. The loop free energy is thus given 
by / oc clnn. Within the framework of the Fokker-Planck equation, the probability 
distribution of finding a loop of size n at time t, P{n, t), is given by 

''^"•'^-ol\z^ir]nn.t), (2) 



dt dn ^ 

where D is the diffusion constant. Here we have taken the continuum limit and assumed 
the dynamics to be over-damped. This equation has to be solved with the boundary 
condition P(0, t) = and initial condition P(ra, 0) = 5{n — nQ). The survival probability 
of the loop is then given by G{nQ., t) = J dnP{n, t). 

Within the scaling approach the survival probability is written in the form 

G{no,t)=g{Dt/n^,) , (3) 

with z = 2. In Appendix 1 we show that the asymptotic behavior of the scaling function 
for small and large values of the argument is 

g{x) ~ 1 for x < 1 (4) 

g{x) ~ X ~ for X ^ 1 . (5) 

The autocorrelation function (Eq. ([T])) may thus be written as 

ff nl~^g{Dt/nl) duo 

r^ i-c 7 
Ji ^0 duQ 



Git) ^ '' '-o^^v^-/-o;--^o ^g^ 
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where the system size A^ is taken to infinity in the thermodynamic hmit. We first 
consider the long time behavior for c < 2. In this case the integrals are controlled by 
the upper limit A^, where g{Dt/n^) ~ 1. Both numerator and denominator diverge as 
N'^~'^ so that C{t) ~ 1 for t ^ 1. On the other hand for c > 2 both integrals are 
independent of the upper limit. Changing variables to y = riQ/yDt yields 

C{t) ^ [Dtf-^l'^- / y'~^g{y-')dy , (7) 

where (rio) is the average loop size. The asymptotic behavior of g{y~^) at small y (Eq. 
([5])) implies that the integral converges for t -^ cxd, yielding C{t) ~ t^'^^"^. Hence 

This expression suggests that measuring C(t) at criticality may be used to determine 
the entropy exponent c. In particular it can be used to distinguish between the case 
of a continuous transition (c < 2), where C{t) = 1, and a first order phase transition 
(c > 2), where C{t) decays to zero at long times. 

In the above analysis it is assumed that the loop is at local equilibrium at any 
given time. For this assumption to be valid, the survival time of large loops has to be 
much longer than its equilibration time. A typical survival time of a loop of length n 
scales as n^. On the other hand, simple models for the dynamics of microscopic loop 
configurations, which are usually based on diffusion processes, yield relaxation times 
which also scale as n^. Thus the two typical times scale in the same way with the 
loop size, and it is not a priori clear that during the evolution of the loop it is at local 
equilibrium. Note that off criticality the loop size changes linearly in time and therefore 
the assumption of local equilibrium is clearly not valid. 

In the following we introduce and study a model for the loop dynamics. This is 
done in two steps: First, we consider the simpler case of c = 3/2 discussed in [17]. We 
then generalize this approach to arbitrary values of c. We find strong evidence that the 
local equilibrium assumption holds asymptotically for the model. It is thus argued that 
within the model the local equilibrium assumption is valid. 

2.2. Microscopic Dynamical Model for c = 3/2 

In this section we introduce and analyze a simple model of loop dynamics corresponding 
to c = 3/2. Within the model, the loop is described by a fiuctuating interface (or a 
string), interacting with an attractive substrate in d = 1 + 1 dimensions. Here the 
interface height variable corresponds to the distance between complementary bases. The 
interface configurations are those of a restricted solid on solid (RSOS) model defined as 
follows (see Fig. ([T])): Let hi = 0, 1, 2 . . . be the height of the interface at site i. The 
heights satisfy \hi — /ij+i| = ±1. Consider a loop between sites and n (where n is 
even) as shown in Fig. ([1]). Outside the loop the interface is bound to the substrate so 
that h_2k = hn+2k = and h-2k-i = hn+2k+i = 1 for A; = 0, 1, . . .. Inside the loop the 
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heights hi...hn-i can take any non-negative value which is consistent with the RSOS 
conditions. For simphcity we allow only one end of the loop to fluctuate while the other 
is held flxed. This should not modify any of our results, since the dynamics of the two 
ends of long loops are uncorrelated with each other. We consider a random sequential 
dynamics in which the loop conflguration and its length are free to fluctuate. Thus the 
dynamical moves are as follows: 

hi ^ hi ± 2 with rate 1 for sites 1 < i < n — 1, (9) 

as long as the resulting heights are non-negative and the RSOS condition is satisfled. 
For i = n the loop length is changed according to the rules 

n —^ n + 2 with rate a/4 

n ^ n — 2 with rate a , (10) 

where n can decrease only if hn-2 = 0. At the other end the height is flxed, ho = 0. It 
is straightforward to verify that the number of conflgurations of a loop of size n is given 
by 2"'/n'^ with c = 3/2 for large n. This is a result of the fact that the number of walks 
of length n in d = 1 + 1 dimensions is 2" and the probability of flrst return is n~^/^. The 
ratio, 1/4, between the two length changing processes in Eq. flTOj) is chosen such that 
in the large n limit the loop is not biased to either grow or shrink. This corresponds to 
the model being at the denaturation transition point, which is determined by equating 
the free energies of the pinned segment and that of the open loop. Combining this with 
detailed balance yields the ratio between the rates. The parameter a determines the 
rate of the length changing processes: a = corresponds to the dynamics of a loop of 
a flxed length. As a is increased the length changing processes become faster. In the 
following subsection this model is generalized to include a power law potential between 
the interface and the substrate. This will allow us to study other values of c. 

In a realization of this dynamics one of the n + 1 attempts deflned above, Eqs. ([9]) 
and (fTOl) . is chosen at any given time. Of these, n — 1 are attempts to update the height 
at sites 1, 2, . . . , n — 1. The other two are attempts to update the position of the edge 
by a move either to the right or to the left. One attempted move of the edge defines a 
Monte Carlo sweep. 




Figure 1. A typical microscopic configuration of tlic loop in the RSOS model. Dashed 
lines indicate possible dynamical moves of the interface. 

In order to test the validity of Eq. ([2]) we compare its predictions with results 
obtained from numerical simulations of the model above. In the numerical simulation 
we find good data collapse, when plotted against t/ng with z > 2, depending on the 
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value of a, rather than the expected Fokker-Planck value z = 2. With these modified z 
exponents the survival probability agrees well with the results obtained from the discrete 
version of the Fokker-Planck equation. The results are summarized in Fig. ([2]) where 
the survival probability is plotted as a function of the scaling variable t/n^'^ and tjri^^'^ 
for a = 1 and a = 0.1 respectively, for several values of the loop size no. The question 
is whether the discrepancy in the value of z is a result of a finite size effect or does it 
persist in the large n^ limit. For the Fokker-Planck equation to properly describe the 
system it is essential to show that z approaches 2 in the large tlq and t limit. 
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Figure 2. Data collapse of the survival probability (averaged over 4 • 10^ realizations) 
for some values of ng with (a) a = 1 and z — 2.2 , and (b) a — 0.1 and z — 2.07. The 
line corresponds to a numerical solution of Eq. ^ . 



In the following we argue that in fact the value z = 2.2 in the case of a = 1 (and 
z = 2.07 for a = 0.1) is a result of finite size effects. For large systems the value z = 2 
is expected to be recovered. To check this point we calculate numerically the variance 
of the loop size 

w\t) = {{n{t) - {n{t))r) , (11) 

where (■) denotes an average over realizations of the dynamics. In order to evaluate 
the temporal growth of w'^{t) we define a variable cr^{t) which takes the value 1 if the 
length of the loop increases at time t and otherwise. Similarly, we define o"_(t) and 
cro(t) for steps which decrease the loop size and steps in which the loop size does not 
change, respectively. Clearly (J+{t) + cr_(t) -|- cro(^) = 1- The dynamics of the chain, Eq. 
(ITOl) . implies that in the limit of large uq one has 

(or+(t)) = ((T_(t)) = a/8 ■ (ao(t)) = 1 - «/4 , (12) 

where a = a/max{l,a} in accordance with the random sequential dynamics. Denoting 

U{t) = <J^{t) — <J-(t), it is easy to see that 

Aw\t) 



At 



w\t)-w\t- 1) 
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t-i 



4(f/(tf)+85^(f/(r)f/(t)) 



(13) 



T = l 



where 

(f/(r)t/(t)) = K(rK(t)) + (a_(rK(t)) 

-(a_(rK(t))-K(rV_(t)). (14) 

It is evident that a loop increasing step at time t, (cr_|_(t) = 1), is uncorrelated with 
steps which took place at time r < t. Thus (cr+(r)cr+(t)) = (o"_(r)(T+(t)) = a^/64. 
Numerically we find (o"_(r)o"_(t)) = a^/64 (see Fig. ([3])). Using these results we finally 
obtain 

Aw\t) 



At 



a 



t-i 

r=l 



[(^+(T)^-(t))c] 



(15) 



with (a+(r)a_(t))c = (a+(r)a„(i)) — a;^/64. Numerical simulations of the dynamics 
show strong correlation between a_|_ (r) and cr_ (t) with an algebraic decay in t — r (see 
Fig. (I3l)). It is interesting to note that the dynamics of the chain induces such long 
range temporal correlations between steps of the edge mediated by the loop dynamics. 
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Figure 3. Correlation functions of the a variables as obtained by averaging over 
1.9 • 10^ realizations, for no=4000. 



By extrapolating the sum on the right hand side of Eq. ( 1151) using the asymptotic 
form B{t — t)~"' with B ^ 0.015 and 7 ~ 1.2, deduced from Fig. ([3]), we find that 
the sum converges to a non-zero value. This is demonstrated in Fig. (j4]) for a = 1 
and a = 0.1. For example, in the case a = 1 the sum converges to ~ 0.84 < a = 1 
indicating that w'^{t) ^ 0.16t at large t, which in turn yields z = 2. The slow power-law 
convergence towards the asymptotic value implies that it may require large systems to 
observe the long time behavior of Eq. ([8]). 
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Figure 4. The diffusion coefficient D = Aw^(i)/Ai , as calculated by Eq. (jlSp . 
averaged over 270, 000 runs with np = 2000 for (a) a = 1 and (b) a = 0.1. The slow 
decay of D can be easily observed. 



2. 3. Microscopic Dynamical Model for Arbitrary c 

In this section we generalize the model of the previous section to consider the case of 
arbitrary c. This can be done within a. d = 1 + 1 dimensional model by introducing a 
repulsive interaction between the substrate and the interface. Taking an interaction of 
the form A/h?, where h is the distance between the interface and the substrate and A 
is a constant, results in an equilibrium weight of a loop of the form 1/n'^. The exponent 
c is related to the interaction strength A [221 123] . 

In order to derive the relation between A and c one notes that at the critical 
point the distribution of the distance between the interface and the substrate decays 
algebraically at large distances, Q{h) ~ 1/h'^. It has been shown that for an interface 
model for which self-avoiding interactions play no role, the exponent k is related to the 
loop exponent c by [21] 

c=(k + 3)/2. (16) 

We proceed by introducing a specific model and evaluate k in terms of the interaction 
parameter A. One then obtains the loop exponent c from Eq. flT6|) . We consider an 
RSOS interface model with the Hamiltonian 

A, 



H{hi,h2,...,hn) = ^ 



-£Sh,fi + 



w 



s, 



hi,0) 



:i7) 



where as before, hi = 0, 1, 2..., and hi — /ij+i = ±1. In this Hamiltonian e > represents 
the binding energy between the substrate and the interface. To evaluate Q{h) we write 
down the eigenvalue equation of the transfer matrix corresponding to the Hamiltonian 
Eq. (fTTI) . For h > 1 the equation is 

e-PA/h^Mj^_^ + e-^^/''^,+, = AvI/, . (18) 
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Here A is the eigenvalue and ^^ are the components of the eigenvector. The distance 
distribution is given by Q(/i) oc \1/^. At criticahty the eigenvector component, at large 
/i, has a form \&(/i) = l/h'^/'^. By using this form in Eq. ( ITSj) we find the relation 

[5A=-k{k + 2) . (19) 

8 

Combining this with Eq. (fT6l) yields 

/3v4= -(2c-3)(2c-l) . (20) 

8 

We now use the model, Eq. (fT7|l . to study numerically the dynamics of a loop with 
c 7^ 3/2. The dynamics of the model is similar to that introduced in Sec. 2.2, but with 
the transition rates of Eq. ([9]) modified according to the Hamiltonian Eq. (IT7I) . Namely, 
the updating rates are given by 

hi —^ hi + 2 with rate 1 

h,^hi-2 with rate e-^^l^^^"^)-^-^-^) ^ ^21) 

as long as the resulting heights are non-negative and the RSOS condition is satisfied. 
The presence of the long-range interactions also changes the ratio between the rates by 
which the loop grows {R{n — » n + 2)) and shrinks {R{n + 2 ^ n)). This ratio is given 
by 



R{n^n + 2) _ ., 



(22) 



R{n + 2^n) 

The critical temperature is found by equating the free energy of the loop with that of 
the bound segment. This yields 

4 = g-/3(A-.) ^ (23) 

where A — e is the energy of a pair of sites in the bound segment and 4 is the statistical 
weight of a pair of sites in the open loop. Combining this with Eq. fl22l) gives 

n — i> n + 2 with rate ae~^'^/4 

n —>■ n — 2 with rate a . (24) 

We have simulated the dynamics of Eqs. fl2Tl) and fl2^ for A = 0.15 and A = 0.5. 
These values of A correspond to c ~ 1.74 (< 2) and c ~ 2.12 (> 2) respectively. In Fig. 
([5]) we present the loop size distribution for these two values of the parameter A. The 
resulting c values fit well with the predictions. 

In studying the survival probability of a loop we follow the same approach which 
was applied in the previous section for c = 3/2. Similar results were obtained for the 
case of c > 3/2. In Fig. ([H]) we present the survival probability as obtained from 
numerical simulations of the model. We find good data collapse, but again with a 
modified exponent z = 2.2 for a = 1. The scaling function fits well with that obtained 
from a numerical integration of a discrete version of Eq. ([2]). 

We have also calculated the step-step autocorrelation function as for the case 
c = 3/2 and found similar results. In particular, we find that the exponent 7 seems to 
have a weak dependence on A, with 7 ~ 1.4 for both A = 0.15 and A = 0.5 (Figures 
not shown). 
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Figure 5. The loop size distribution for A — 0.15 and A — 0.5 as measured in 
numerical simulations, and theoretical resulting exponents c ~ 1.74 and c ss 2.12 
respectively. The theoretical curves show good fit to the measured data. 
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Figure 6. Data collapse of the survival probability (averaged over 4 • 10^ realizations) 
for some values of n^ with z — 2.2, with (a) c = 1.74 {A = 0.15) and (b) c = 2.12 
{A — 0.5) . The line corresponds a numerical solution of a discrete version of Eq. ([2]) 
with corresponding values of c 



3. Many loops model 

In the previous section we analyzed the dynamics of a single loop. We found that 
it is well described by the Fokker-Planck Eq. ([2]) for asymptotically large loops. In 
the present section we extend this model to consider interaction between loops. This 
is done by considering a chain composed of an alternating series of loops and bound 
segments. Each loop and bound segment is characterized only by their respective length. 
In contrast to the study of the dynamics of a single loop here no internal degrees 
of freedom are associated with a loop. Within this model loops evolve by growing, 
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shrinking, splitting, merging, together with creation and annihilation processes. The 
rates of the various processes are chosen so that the system evolves to the equilibrium 
loop length distribution at large times. While the choice of rates is not unique they are 
taken to be compatible with the single loop dynamics whenever applicable. A similar 
approach has recently been applied to study dynamical features such as the approach 
to equilibrium near the denaturation transition [18j. From this analysis we extract the 
behavior of the autocorrelation function of a base-pair inside a dsDNA where many 
interacting loops coexist. 

3.1. Definition of the Model 

The DNA configurations can be represented by an alternating sequence of bound base- 
pairs and loops. We denote by [k] a bound segment with length k and (/) a loop 
of length /, with k,l > 0. A given configuration of the DNA is thus represented by 
[^i](^i)[^2](^2) • • •• In terms of these variables the dynamics of the model is defined by 
the following rates: 

• Motion of a loop edge. This corresponds to the same processes which were 
considered in the dynamics of an isolated loop in the previous section. 

[A;](/) ^ [A; -!](/ + 1) with rate (^)' 

[k - l]{l + 1) ^ [k]{l) with rate 1 

(/)[A;] ^(/ + 1)[A;-1] with rate (^)' ^ ^ 

(/ + 1) [A; - 1] ^ (/) [k] with rate 1 

These processes are executed as long as the lengths of the resulting loops and bound 
segments are non-zero. 

• Splitting and merging of loops 



(/i + /2 + l)^(/i)[l](/2) with rate ^ (^ 



'i'2 J (26) 
(/i)[l](/2) ^ (/1 + /2 + I) with rate 1 

In addition we consider creation and annihilation of loops. 

• Creation and annihilation of loops 

[fci + fc2 + l]->N(l)N with rate ^(^y^^ .^7) 

[A;i](l)[A;2] -^ [A;i + A;2 + 1] with rate 1 

Here ctq is the cooperativity parameter, and C(c) = Xl-^i "'~'^- I^ is straightforward to 
verify that the choice of rates satisfies detailed balance with respect to the equilibrium 
weight for the loop sizes at criticality P{n) = o'ottt- 

3.2. Numerical Simulation 

To check that indeed interactions between loops do not modify the asymptotic behavior 
of the autocorrelation function C(t) we simulate the model Eq. (I25l) - (p7|) . We use 
the experimental relevant value uo = 10~^ [19] and consider a DNA length of 100, 000 
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base-pairs. The autocorrelation is evaluated by monitoring the state of 1000 base-pairs 
uniformly distributed within the DNA. Fig. ([7]) shows the results for various values of 
c along with the theoretically expected slopes. 
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Figure 7. Normalized autocorrelation functions as measured in the simulation of the 
many loops model, for a^ = 0.0001, L = 100000 and 50000 repetitions. The thin lines 
indicate the expected behavior of the autocorrelations for the appropriate values of c. 



While the results for large values of c agree well with the theory, there is a systematic 
deviation from the predicted slopes for smaller values of c close to 2. These deviations 
could be attributed to the finite length of the simulated system. For example it is clear 
that for c < 2 the autocorrelation function of a finite system decays to zero at long 
times rather than remaining constant. This is due to the fact that there is an upper 
cutoff on the loop size available. Only for an infinite system C{t) is expected to remain 
constant (= 1) at long times. In order to check this point we introduce an upper cutoff 
Nmax to the loop size in the equation for the autocorrelation function 



^^^, 



Cit) 



EZ:: Pe,ino)noGin,,t) 



^M 



El'lT Pe,ino)no 



(28) 



The loop size N^ax is chosen so that it appears roughly one time during a run. For runs 
which are not too long this can be estimated using a^LRP^Nmax) = 1, where L is the 
system size and R is the number of Monte-Carlo repetitions performed. In Fig. ([H]) the 
results of the simulations are compared with the theoretical expression, Eq. fl28|) . which 
is summed numerically. 

In summary, we find that our scaling predictions are generally confirmed by the 
numerical simulations of the many loops model. However, for values of c close to 2, 
deviations are found. These seems to be related to finite size effects. 



Dynamics of DNA Melting 



14 




Figure 8. Same as Fig. ([5]), witii tiic numerical calculations of the sum with a cutoff. 



4. Conclusions 



In this paper the dynamics of loops at the denaturation transition was studied both 
within a single loop model and a many loop approach. In particular, special care was 
given to the applicability of the Fokker-Planck equation. It was shown that the long- 
time decay of the autocorrelation function of the state of complementary bases (closed 
or open) is sensitive to the value of the loop exponent. In particular, for c < 2 it remains 
finite while for c > 2 it decays as t^~'^''^. 

Throughout the paper we have considered homopolymers where the binding energy 
between different base-pair is constant. In typical DNA molecules the binding energy is 
not homogeneous. While a preliminary treatment of the effects of disorder was given in 
[TT] . it remains an important and interesting question. 

Acknoledments: The support of the Israeli Science Foundation (ISF) and the 
Albert Einstein Minerva Center for Theoretical Physics is gratefully acknowledged. YK 
also acknowledges support by the US-Israel Binational Science Foundation (BSF). 



5. Appendix I: Asymptotic behavior of the return probability 



In this appendix we derive the asymptotic behavior of the survival probability 
corresponding to the Fokker-Planck equation 



dP{n,t) _^d 
dt ^ dn 


'c d' 
n dn 


P{n,t) 






(29) 


with the boundary conditions 








P(0,t) = ; P 


(oo,t) = 


; P(n,0) = 


= 5{n- 


-no) . 


(30) 
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To do so, we first perform a Laplace Transform 



oo 

st \ 



P{n,s) = / e-''P{n,t)dt (31) 

Jo 

to obtain 

sP{n, s) - 5in - no) =D^\^ + ^-1 P{n, s) . (32) 



Integrating over a small interval around no yields 

dnP<(n)\ - dnP>(n)\ =— (33) 

" ^ V / 1 n=no ^ ^ ' I n=no JJ ^ ' 

where P<,{n) and P>{n) are the solutions of Eq. fl32l) for n < uq and n > uq respectively. 
By defining x = ^ sjDn and P(n, s) = {Ds)~'i f{^/s/Dn) Eq. ([32D becomes: 



r(x) + -fix) - (i + 4) /(^) = (34) 

X \ x^ / 

i solution 

fix) = Ax^Ii+c (x) + Bx^Ki+c (x) . (35) 



X 

which has the solution 

1 — c , , 1 — c 

i±c(x) + Bx 2 _ 

2 2 

Here Jj, and i^;^ are modified Bessel functions of the first and second kind 
Using their asymptotic behavior and the boundary conditions (130|) we find B = 
for X < Xq = ^ sjDn and A = for a; > x^. Denoting x< = min{x,xo) and 
a;> = max{x, Xq) and using Eq. fl33l) the Laplace Transform of the loop size distribution 
is given by 

/ \ (l-c)/2 
- [fj I(l+c)/2{x<)K(i+,),2{x>) 

P{s,x) = ^^ ^^ ^ (36) 

^D^ \ I[l+c)/2i^o)K^l+c)/2{xo) - /(i+c)/2(a;o)i^(i+c)/2(3^o) 



Using standard methods [2Z1 we integrate this expression to find the Laplace transform 
of the survival probability 



G{s,xo) 



POO 

/ P{s,x)^/D/s dx 
Jo 



c-l 



Ki+c{xo) [ I^xq) - p.'.!^" ^l^ ) + Ii±^{xq)Ki^{xq) 



r((i+c)/2) 

(37) 



s {ri+c{xQ)Ki+c{xQ) - Ii+^{xq)K'i+,{xq) 

V 2 2 2 2 

The asymptotic behavior of G{s^nQ = ^D/ sxq) for long and short times can be 
extracted from the behavior of the Bessel functions. For small s Eq. (1371) turns into 



c-l 



—, , $(c)s^" c < 1 , , 

^t--).!^^,,^,^^^;, (38) 

where $(c) is a constant which depends on c. From this we can extract the asymptotic 
form of the survival probability for long times: G{no,t ^ n^/D) = g{^ = Dt/nl ^ 

1 + c 

1) ~ ^~^~ ■ The behavior for short times can be obtained in a similar fashion, yieding 
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g{^ <^ 1) K, 1. In sum, we find that the survival probabihty for a loop of initial size tiq 
has the scaling form 

G{n,,t)=g(^\ , (39) 

V^o / 
With the asymptotic behavior 

g{^ > 1) ~ r^ 

g{^ « 1) ~ 1 . (40) 
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